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Abstract 

Conditions for almost sure extinction are studied in discrete time 
branching processes with an infinite number of types. It is not assumed 
that the expected number of children is a bounded function of the par- 
ent's type. There might also be no integer m such that there is a lower 
positive bound, uniform over the ancestor's type, for the probability that 
a population is extinct at the m-th generation. A weaker condition than 
the existence of such an m is seen to lead to extinction almost surely if 
the sequence of expected generation sizes does not tend to infinity. Some 
criteria for a positive probability of nonextinction are given. Examples 
are provided by extending to our setting two applications, namely Leslie 
population dynamics and processes arising in continuum percolation in 
which the offsprings follow Poisson point distributions. 
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1 Introduction 

We study the subject of giving conditions for the extinction probability to be 1 
in discrete time branching processes with an infinite number of types. Known 
results suppose that the expected number of children is a bounded function of 
the type of their parent (Harris, 1963; Mode, 1971). Here we do not prevent 
types from being arbitrarily large. This means that there might be no uniform 
bound over the types for the expected number of children, nor a number m of 
generations such that there is a uniform positive lower bound for the probability 
that an individual will have no descent at the m-th generation. A weaker con- 
dition than the latter is introduced in Section 2 and it is shown to be sufficient 
for almost sure extinction if the sequence of expected generation sizes does not 
tend to infinity, even if these expectations fluctuate without an upper bound. 
Existence of a strictly dominant eigenvalue of the integral operator associated 
to the offspring means is not assumed in our propositions. Focusing on the 
expected number of individuals in subsets of types that are in some sense larger 
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than the type of the initial individual, criteria for the extinction probability to 
lie below 1 are obtained in Section 3. Section 4 provides examples by extending 
to unbounded types a couple of well known applications. The possibility of con- 
structing counterexamples shows that the sufficient condition for almost sure 
extinction in Section 2 is not always necessary and leaves a criticality problem 
open. 

2 Almost sure extinction 

Let Zq,Z\,... be a discrete-time branching process with types in a space X, 
countable or not. In the uncountable case, X is an Euclidean space and Z , Z\, ... 
a sequence of point distributions that follows the branching process definition 
as given in Harris (1963), Chapter 3. The number of individuals with type 
in a measurable set A C X at the n-th generation will be denoted by Z n (A). 
The symbols P X Q, E x and q x will stand for probabilities, expectations and the 
extinction probability if the process is initiated with a single individual of type 
x 6 X. 

Condition 1 For any number of children k > 0, there exist a number of 
generations m(k) and a real number q(k) > such that 

mi P x (Z mW (X) = 0|1 < Z X {X) <k)> q(k). 

xEX 

Proposition 2 If Condition 1 holds and liminf„ E x Z n (X) < oo, then 
q x = 1. 

Proof. By Fatou's Lemma, liminf n E x Z n (X) < oo implies 
E x liminf„ Z n {X) < oo, i.e. J2k>o -fx(hminf„ Z n (X) > k) < oo. Thus given 
e > 0, there is a positive integer k € such that Pa,(liminf„ Z n (X) > k e ) < e. 

So P x (Z n (X) > k e for all but finite n's) < e and therefore 

P x (Z n (X) < k e ,i.o.) > l-e. 

We shall see that for any positive integer k, P x (l < Z n (X) < k,i.o.) — 0. This 
will complete the proof because then for all e, P x (Z n (X) = 0, i.o.) > 1 — e and 
the extinction probability must equal 1. Let 

N {k) = 

and for i = 1,2, 

Ni(k) = min{n > N t ^{k) - 1 + m(k) : 1 < Z n (X) < k} 
with min0 = oo. Now, 

P X {1 < Z n (X) < k,i.o.) = P x (f) t {N t (k) < oo}) = lim P x (Ni(k) < oo), 
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so it is enough to prove that P x (Ni(k) < oo) < (1 — q{k) k ) 1 . 
For i = we have 

P x (N (k) < oo) = 1 = (1 - q{k) k f 

and for i > 

P x (N i+1 (k) < oo) < P x (Ni(k) < oo,Z N . (k) _ 1+m(k) (X) > 0) 

OO 

= Px(Z n -i+m(k)(X) > 0\Ni(k) = n)P x (N t (k) = n). 

n=l 

But given Ni(k) = n, there are at most k individuals at generation n. These 
must be the children of at most k individuals out of generation n — 1, which have 
at most k children each. Hence by Condition 1,1 — q(k) k is an upper bound for 
P x {Z n -i +m (k)(X) > 0\Ni(k) = n) yielding 

OO 

P x (N l+1 (k) < oo) < (1 - q(k) k ) Px{Ni(k) = n) 

71=1 

= (l- g (fc) fe )P x (7V i (fc)<oo). 
The desired result follows by induction. 

3 Positive probability of nonextinction 

Let Zq, Z\, ... be a discrete-time branching process with types in a countable or 
uncountable set X as in the preceding section. Given a type x and a positive 
integer n, by the n-th generation descendants of x we will mean the individuals 
in the n-th generation of a population that begins with a single individual of 
type x. Let us consider the following order relation in the set of types. Given two 
types x and y we shall say that x is larger than y if for every measurable subset 
A of types and any n (n = 1,2, ...), the number of n-th generation descendants 
with types in A is stochastically larger for x than for y. 

Proposition 3 Let y be a type in X. Suppose that there exist a measur- 
able set A whose elements are larger than y and a positive integer i such that 
EyZi(A) > 1. Then for any x larger than y we have q x < 1. 

Proof. Let So, Si, ... be a single-type branching process defined by 

S = l 

and for k = 0, 1, 

P(S 1 =k)=Py(Z i (A)=k). 

As EyZi(A) > 1, this process is supercritical. Now consider Z ,Zi,Z 2 i, the 
multitype process at those generations that are multiples of i, beginning with a 
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single individual whose type is x. Since x and all elements of A are larger than 
y, Z ni (A) is stochastically larger than S n for any n. As Z ni {X) > Z ni (A), we 
have q x < 1. 

Corollary 4 If for a type x there exists i such that E x Zi({x}) > 1 then 
q x < 1. 

Corollary 5 Suppose that there exist y, i and A as in the proposition. If 
there is a type x, an integer j and a set B C X whose elements are larger than 
y, such that E x Zj(B) > 0, then q x < 1. 

Corollary 6 If there is a smallest type s in the sense that every type x is 
larger than s, and E s Zi(X) > 1 for some i, then ^ < 1 for every x. 

4 Examples 

Before going into the examples it should be remarked that in all the processes 
we are dealing with, the expectations E x Z n (A) verify 



for n = 1,2,..., as proved in Harris (1963), Chapter 3, and in Mode (1971), 
Chapter 6. The set A is any measurable subset of X and dE x Z\{y) means that 
the integration is performed with respect to the measure defined by E X Z\{S) 
for any measurable subset S of X. Under discrete interpretation, the right 
hand side of the equality becomes the result of rising the first generation mean 
matrix to the (n + l)-th power and taking the sum of the elements in the row 
of x that lie in the columns that correspond to types in A. In the continuous 
case, suppose that k(x,y) is a nucleus such that for any measurable A C X we 
have E X Z\{A) = J A k(x,y)dfj,(y) for some measure fi on X. Then the formula 
yields the iteration of the integral operator K associated to k and fj,, allowing 
us to write E x Z n+1 (A) = j ' A K n {k{x,y))d^{y). 

In all our examples we shall both have 



for any m. That the expectations are not bounded will be evident. In order to 
prove the second property, the following proposition will be useful. 

Proposition 7 Suppose that for any positive integer h, 
mi x( zx Px{Zi(X) < h) = and that there exists a > such that for every 




sup E X Z\{X) = oo 

xEX 



as 



inf P x (Z m (X) = 0) = 
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x G X, P X {Z X {X) > 0) > a. Then mf xeX P x (Z m (X) = 0) = for any number 
of generations m. 

Proof. Let {x h : ft = 1, 2, ...} C X be such that P Xfc (Zi(X) < ft) -> when 
ft — > oo. For any m, 

P Xh (Z m (X) = 0) 
< P Ih (Zi(X) < ft) +P Xh (Z m pO = 0|Zi(X) > K)P Xh (Z x (X) > ft) 
<P Ch (Z 1 (X)<ft) + (l- a m - 1 ) ft . 
So liminf/j P Xh {Z m {X) = 0) < and the proposition is proved. 

4.1 Leslie dynamics 

We shall consider the infinite types analog of the finite-types branching process 
studied by Pollard (1966), which is the stochastic version of the deterministic 
population dynamics described by Leslie (1945). The set of types will be here 
{xo, xi, ...}. If we adopt the interpretation of types as age classes, then the 
individuals of type xo are the youngest ones. An xo may bear individuals of 
type xo or may produce an individual of type x\ by surviving enough time to 
enter in the next age class. In general Xi (i = 0, 1, ...) may only produce Xo's 
and eventually an Xj+i, so that the infinite mean matrix has only two positive 
entries at each row, E Xi Zi({xo}) and E Xi Z 1 ({xi + i}). For simplicity, both the 
individuals x as the x i+ i arising from an Xj will be called its children. 

Our examples will have first generation means E Xi Z\({xo)) = p 1 ^ 1 and 
E Xi Zi({xi + i}) = p 2 for p in (0, 1), so the infinite mean matrix is 

I p p 2 ... \ 

I p 2 

p- 1 p 2 

p- 2 p 2 

\ - - / 

By what we noticed at the beginning of this section, the first row of the n- 
th power of this matrix consists of the n-th generation means E Xo Z n ({xo}), 
E Xo Z n ({xi}) It can be seen by induction over n that these entries are 

f (I/2)(2p)"(p/2) 1 if 0<i<n 
E X0 Z n ({ Xi }) = I p 2n if i = n . 

{ if i> n 

In particular, E Xo Z n ({x }) = (l/2)(2p) n so by Corollary 4, q XQ < 1 if p > 1/2. 
By a simple case of Corollary 5, q Xi < 1 for any i. 

In case p < 1/2 the mean number of individuals of any type at the n-th 
generation (n = 1, 2, ...) in a population that started with a single individual of 
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type xq is bounded by a constant. Indeed 



E X0 Z n (X) =Y,E X0 Z n {{ Xi }) = (l/2)(2p)" ^ Pl ' +p 2n < 1. 

Thus liminf„ £7 Xo Z„(X) < oo is satisfied. But if p < 1/2 we can bound the 
expected total number of descendants in all generations because 

OO _| OO OO 

E xo z n (x) < — x> P )" + E p 2n < °°- 

71=1 n=l n— 1 

So £ Xo X^=i Z n (X) < oo and therefore = 1. 

We shall now consider two different offspring distributions behind our means. 
In both examples, the number of a^-type children of an Xi parent (i = 0, 1, ...) 
will follow a Poisson distribution with mean p x ~ % and the number of Xi+i- 
type children will be 1 with probability p 2 and with probability 1 — p 2 . The 
difference will lie in the joint distribution for these marginals, as described below. 
It can already be seen that ini x£ x P x (Z m (X) = 0) = for any m because having 
these marginal distributions implies the hypotheses of Proposition 7. 

The first choice of distributions for the children of the a^'s (i = 0, 1, ...) is 
one that satisfies Condition 1. We take for each Xi the joint distribution that 
makes the number of children of type x independent of the number of x i+ i 
children. This family satisfies Condition 1 for m(k) = 2 because 

P Xi (Z 2 (X)=0\l<Z 1 (X)<k)> 

> P Xi (Z 2 (X) = 0|Zi({a: i+ i}) = 0, 1 < Z X (X) < k)- 
■P Xi (Z 1 {{x i+1 }) = 0|1 < Zx[X) < k) 

where 

p Xi (z 2 (x) = olZi^i+i}) = o, i < z x {x) <k)> (P^izAx) = o)) fe , 

while for any integer h between 1 and k, 

P Xi {Z 1 {{x i+1 }) = Q\Z 1 {X) = h) = 

(l-p^e-^y-y/fe! _L_ 

~ (l-p 2 )e-P 1 - l (p 1 -') h /h\+p 2 e-P 1 - l (p 1 -*) h - 1 /(h-i)\ > 1 + jfe" 

Thus 

inf ■ P Xi (Z 2 (X) = 0|1 < Z,{X) <k)> {P XQ {Zi{X) = 0)) k -^r. 

Since it has been shown that liniinf n E Xo Z n (X) < oo for p < 1/2, we have 
q Xo = 1 by Proposition 2. 
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The second family of distributions does not satisfy Condition 1 because there 
is enough positive correlation between bearing few children and having one very 
prolific among them. A way to achieve this is by requiring that at least for 
k = 1 and for infinite i's, 

P Xi (Z 1 ({x i+1 }) = 1|1 < Z X (X) < k) = 1. 

This is possible because due to the marginal distributions P Xi (Z\(X) < 1) — ► 
when i — > oo. So for any i greater than some i, P Xi (Z\ (X) < 1) < p 2 . Since p 2 is 
the probability of having one Xj+i, we can define for those i's a joint distribution 
that loads all of P Xi (Z\(X) < 1) only on the probability of having no xq and 
one Xj+i. So for i > i, P Xi (Zi({xi+±}) — l\Z\(X) = 1) = 1, as we required. 
Now we are satisfying an analogous hypothesis as that of Proposition 7, namely 
that for any integer h, m{ Xi ^x P Xi {Z%(X) < h\Z\(X) = 1) = 0. By this and 
the fact that there exists a > such that for every Xi, P Xi (Z\(X) > 0) > a, an 
analogous proof as that of Proposition 7 gives 

inf P Xi (Z m (X) = 0\Z 1 (X) = l) = 

XiGX 

for any number of generations m. This means that Condition 1 is not satisfied 
because it fails for k — 1. Since it was already proved that q Xo = 1 for p < 1/2, 
a consequence of this example is that the hypotheses of Proposition 2 are not 
necessary for the extinction probability to be 1. For p = 1/2, the criticality 
problem remains here unsolved. 

Finally it can be said about the case p < 1/2 for i = 1, 2, that q Xi = 1 if 
q Xo = 1 by types communication. In fact, q Xi < 1 implies q Xi _ 1 < 1. 

4.2 A non recurrent variant 

The following modification of the unbounded-types Leslie scheme provides an 
example of the application of Proposition 2 in a non recurrent case. Let X = 
{xo,Xi, ...}, p e (0, 1) and the infinite mean matrix be 

/ p ... \ 

1 p 

p- 1 p 

p- 2 p 

V - ••• / 

i.e. the only positive expectations are E Xi Zi({x n }) = p 1 ^ 1 for i = 0,1,... 
and E Xi Zi({xi + i}) = p for i = 1,2,.... In a population that begins with a 
single individual of type xi, only types xq and x n+ i can be present at the n-th 
generation and it can easily be seen that 

i^«({xo}) = ]>> 1 

i=0 
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and 

E Xl Z n ({x n+1 }) = p n . 

Hence 

n 

E xl Z n (X) = W<- . 

In the same way as in the preceding examples, we can find a family of 
distributions with the given means that make q Xl = 1 by Proposition 2, and 
another family of distributions can be defined, such that the criticality problem 
is not solved by Proposition 2. Anyway we notice that although the line of 
descent Xi, Xi + i, ... (i > 1) is a source of growing quantities of individuals of 
type x , it is clearly once broken with probability 1. So as q Xo = 1 we have 
q Xi = 1 for any i > 1. 

4.3 Poisson process descent 

Continuum percolation clusters in R d can be dominated by coupling with branch- 
ing processes so that extinction in the branching process implies that the perco- 
lation cluster is finite. This is done by defining a parent-child relation between 
objects that intersect each other in the percolation setting. Extra individuals 
are eventually added to the families so defined, so that the resulting families fol- 
low the independent offspring distributions of a branching process. This method 
appears in Hall (1985) and it is used and referenced in Meester and Roy (1996). 
It is a way to obtain an upper bound for A, the Poisson intensity parameter of 
the objects' centers, such that below this bound the finiteness of the clusters is 
almost sure. 

This application motivates us to analyze the criticality subject in branching 
processes that would dominate continuum percolation clumps of objects whose 
size might be arbitrarily large. We shall consider the case of spheres in R d 
(d > 1) with not uniformly bounded random radii, which will be the types 
in the branching process. Given a sphere in R d , its offspring is determined 
by spheres that have centers in a homogeneous Poisson process in R d with 
intensity A and have iid radii. The children are those spheres that intersect 
the given sphere. The offspring of each of these children is generated by a new 
independent Poisson process and independent radii an so on. It can be checked 
that we both have sup xeX E X Z\(X) = oo as inf xe x Px{Z m {X) = 0) = for any 
m by Proposition 7. 

Let F be the cumulative distribution function of the radii and assume that it 
has an inverse F" 1 : [0, 1) — > [0, oo) when restricted to [0, oo). A way to deduce 
the expression of a mean nucleus is to realize children spheres as homogenous 
Poissonian points with intensity A in R d x [0, 1). The first d dimensions give the 
center of the child sphere and the last one is the image of its radius by F. Let 
u € [0, 1) be a value for the last component of a point in our representation. We 
know that spheres with radius F _1 (u) intersect a parent sphere of radius x if 
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their centers lie inside the parent sphere or outside of it but at an R d distance 
within (x,x + to the parent's center. So in our (d + l)-dimensional 

representation, the points that correspond to F~ 1 (u)-type children of an x lie 
in the <i-dimcnsional sphere that has the same center components as the center of 
the parent sphere, radius x+F^ 1 (u) and last component u. For any x\ <E (0, oo), 
the union of such spheres over u in the interval [F(0),F(xi)] is a solid with 
volume ^ 

/ Xl v d {x + F-\u)) d du= r v d (x + y) d dF(y), 

JF(0) JO 

where Vd is the volume of the unit sphere in R d . Thus we obtained the volume 
of the subset of R d+1 whose points represent those children of an x that have 
last component u G [F(0), F(xi)] or equivalently, type in [0, Xi]. As our Poisson 
process is homogeneous with parameter A, the expected number of children of 
an x with types in [0, Xi] is 

E x Z 1 ([0,x 1 }) = \ f 1 v d (x + y) d dF(y). 
Jo 

Let kd{x, y) = Xvd(x+y) d . If we suppose that the distribution of the radii has 
finite 2d moment (i.e. sphere content has finite variance) then the integral oper- 
ator Kd : L 2 ([0, oo), dF) — > L 2 ([0, oo), dF) defined by kd is compact because kd 
is a Hilbert-Schmidt nucleus. Indeed, / °° k 2 (x,y)dF(x)dF(y) < oo. Since 
Kd is also strictly positive, it has a dominant eigenvalue p > and therefore by 
£7 x Z„([0,oo)) =' f~ ' K2- 1 {k d {x,y))dF{y), the sequence {E x Z n ([0,oo)) : n = 
1, 2, ...} is bounded if and only if p < 1, for any initial type x. 

It will now be seen that for any dimension d and every type x, we have p < 1 
if and only if q x = 1 . To complete the hypotheses of Proposition 2 when p < 1 
it only remains to show that Condition 1 is satisfied. Fixing x > 0, we shall first 
prove that there is a positive bound for the probability of having all children 
with types below x, or more precisely that for any k > 1, 

inf P I (Zi([0,oo)) = ^i([0,x])|l < Zi([0,oo)) < k) > b(k) 

xG[0,oo) 

for some b(k) > 0. By the representation of children as homogeneous Poissonian 
points, we see that given a number of children, each choice of type for a child 
is independent of the others and the probability that the choice falls below x 
when the parent type is x, equals the volume ratio 

fZv d (x + y) d dF(y) = x d J~ x dF(y) + gti(| jjMFjy) 
J °°Vd(x + y) d dF(y) x d + ^ =1 {f )x d ^ / °° tdF(y) 

The integrals are finite by the assumption about the moments, so this probability 
is a continuous and strictly positive function of x. Letting x — > oo gives the 

strictly positive value f£ dF(y). Hence the function has a positive lower bound 
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f3~ and we can take b(k) = Finally, since an individual whose type is in [0, x] 

x X 

has no children with a probability larger or equal than that for x, Condition 1 
is satisfied by 

inf P x (Z 2 ([0,oo)) = 0|1 < Zi([0,oo)) < k) > fiiP-tZ^X) - 0)) k . 

xG[0,oo) x x 

If p > 1, taking x — as the initial type, E o Zi([0, oo)) > 1 for some 
generation i. By Corollary 6, we have < 1 for every x G [0, oo) if we prove 
that is a smallest type. In fact, x > y implies that x is larger than y in the 
sense of Section 3. This assertion can be justified by coupling, constructing a 
process that begins with an individual of type y jointly with one that begins 
with an x. We first assign children to the y by placing points with intensity A in 
the solid of R d+1 that corresponds to y's first generation descent, i.e. the union 
over u G [0, 1) of the R d spheres of radius y + This set is contained in 

the solid corresponding to the first generation descent if the parent's type were 
x because y + F~ x (u) < x + F^ 1 (u). So we can assign to x the same children as 
those of y plus eventually some more in order to follow the children distribution 
of an x parent. In the subsequent generations, x receives the same descent as 
that of y plus the descent of those first generation extra children. This makes 
every n-th generation at any A C X, equal or larger for x than for y. 
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